BROWN-HET-943, |hep-ph/ 9407299 



^1" 

On 



3 
ON 



^ 



X 



High Energy Scattering in 2+1 QCD: A Dipole Picture 



Miao Li0 and Chung-I Tani 

Department of Physics 
Brown University 



rn ' Providence, RI 02912 

> 

A dipole picture of high energy scattering is developed in the 2+1 dimensional QCD, 
following Mueller. A generalized integral equation for the dipole density with a given 
separation and center of mass position is derived, and meson-meson non-forward scattering 
amplitude is therefore calculated. We also calculate the amplitude due to two pomeron 
exchange, and the triple pomeron coupling. We compare the result obtained by this method 
to our previous result based on an effective action approach, and find the two results agree 
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1. Introduction 

One of the most striking aspects of high-energy hadron-hadron scattering is the con- 
tinued increase of the total cross section cry with the energy. Traditional approaches to 
high energy near-forward hadronic collisions invariably involve the notion of a Pomeron 
exchange. There are currently two seemingly conflicting interpretations of Pomeron in 
QCD: One based on perturbative leading log approximation (LLA) |]J and another 
based on nonperturbative (large- iV c and/or phenomenological) considerations J| . Thus 
one of the main puzzles for our understanding of high energy hadronic collisions in QCD 
is the relation, if any, between the perturbative (Lipatov) Pomeron and the nonperturba- 
tive (soft) Pomeron. In a perturbative treatment, a Pomeron loosely corresponds to the 
color-singlet bound state of two (reggeized) gluons. In a nonperturbative treatment, the 
Pomeron is thought to correspond to nonperturbative gluon exchanges having the topology 
of a closed string. It has been emphasized in || that a "two-gluon ladder" also has the 
topological structure of a "cylinder" in the color space. This fact offers the possibility of 
a unified approach to both perturbative and nonperturbative Pomerons in a systematic 
large- N c treatment. 

Two recent studies in QCD at high energies, together with the small- a; results coming 
from HERA, have heightened interests in this area of research. Formally, a high energy 
hadronic collision in the near- forward limit corresponds to the mixing of a "short-distance" 
phenomenon in the longitudinal coordinates with a "long-distance" phenomenon in the 
transverse coordinates. By treating the longitudinal and transverse degrees of freedom 
separately, one could hope to achieve a "dimensional reduction", simplifying QCD at 
high energies to an effective two-dimensional field theory. [[| An interesting attempt in this 
direction has been made recently by Verlinde and Verlinde |7] . The effective theory involves 
several fields and two coupling constants, the original gauge coupling g and an effective 
coupling e 2 ~ g 2 logs. Like previous attempts by Lipatov and collaborators, the emphasis 
is on first considering individual quark-quark scattering through which one reconstructs 
physical hadronic scattering amplitudes. 

The second important development is due to Mueller || 0, who provided an intuitively 
attractive dipole picture for the high energy hard processes. In this approach, which is 
based on a large- A^ c analysis, the basic scattering is between a pair of color dipoles through a 
two-gluon exchange mechanism and the whole complexity of high energy hadronic collisions 
lies in calculating the wave function and the "dipole density" in the transverse space created 
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by a fast moving hadron. Interestingly, the BFKL pomeron is rederived at the stage of 
calculating the dipole density. For related work, see [ffDfl . One possible advantage of this 



new approach is the relative ease in deriving the BFKL pomeron and generalizing it to 
multi-pomerons case. 

In this paper, we carry out a study for QCD in 2+1 dimensions using this dipole 
picture. One of the purposes of this study is to clarify as well as to amplify the large- 
N c cylinder structure, i.e., closed-string like structure, of the Lipatov Pomeron. Since 
much simpler analytic calculations are involved, the resulting physical picture becomes 
clearer. Secondly, since we have previously carried out a similar analysis for 2+1 QCD using 
Verlindes' approach [p+| , this study also helps to check on the validity of this alternative 
approach. One intriguing result of Ref.|ll[] is the fact that the total cross section due 



to one-Pomeron exchange does not grow with energy. We shall reproduce this result 
using this dipole approach and identify its physical origin as due to color-screening in 
2+1 dimensions. Even more importantly, this study potentially could help to clarify and 
contrast the relation of a perturbative approach to that of a nonperturbative soft Pomeron, 
which also relies on a large- iV c expansion. 

In Sec. 2 we first derive an integral equation governing the light-cone meson wave 
function with emission of soft gluons. A physical interpretation for the virtual correction 
factor involved in the integral equation is also provided. We emphasize the surface in- 
terpretation of the dipole picture in Sec. 3, thus providing a possible link between this 
perturbative scheme and other large iV c nonperturbative schemes. The integral equation 
for the dipole density is solved in Sec. 3, and the meson- meson scattering amplitude due to 
one Pomeron exchange is calculated and found, in agreement with results of Ref. ||11|| . We 



proceed in Sec. 4 to derive and then solve the integral equation for dipole pair density. A 
surprisingly simple result is obtained which allows one to calculate the contribution to an 
elastic amplitude due to two Pomeron exchange. The triple Pomeron coupling is discussed 
in Sec. 5. A brief discussion of our results is given in Sec. 6. 

2. The meson wave function 

In this section we shall write down the integral equation governing the generating 
functional of probability amplitudes of soft gluons in a given fast moving meson in 2+1 
dimensions. This is a preparation for deriving an integral equation governing the dipole 



density in a meson state. We adopt the approach of @ and try to stay as close as possible 
with notations introduced in ||. 

Let ij){k, z) be the light-cone wave function of the meson containing no soft gluons, 
where k is the relative transverse momentum between quark and anti-quark and there is 
only one component in our situation, z is the fraction of total longitudinal momentum 
carried by quark. We shall ignore spinor indices, since they are irrelevant in the following 
derivation. The light-cone wave function with one soft gluon is 

k 2 



ip a (k u fc 2 , zi, z 2 ) = 2gT a (ip(k u z±) - ip(fa + k 2 , zi)) 



2.1 
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where k\ is the transverse momentum carried by quark, k% is the transverse momentum 
carried by the soft gluon, z 2 is the fraction of total longitudinal momentum carried by this 
soft gluon and is assumed to be much smaller than z\. The formula is nearly identical to 
that obtained in |J and is schematically represented by diagrams in figure 1. Note that in 
2+1 dimensions, there is only one polarization and we assumed the polarization constant 
be unity in the above formula. A infrared cut-off \i is also introduced. 





Fig.l 
It proves convenient to express the expression (|2.1|) in coordinate space, 

ip a {xQ,x ll x 2 ,zi,z 2 ) = igT a i(j(xo, xi, Z!)[e(x2 - xq) - e(x 2 - x x )], (2.2) 

where e is the step function, i.e., e(x) = 1 if x > and e(x) = — 1 if x < 0. Just as in the 
3+1 dimensional case, the wave function with one soft gluon factorizes: It is the product 
of the wave function of valence quarks and an elementary function. It follows that the 
probability for one soft gluon emission is then 

(2.3) 



g 



'^{xq^x^x^zuz^ = — C F K(x 2 ^x 21 )4> { - \xo,x 1 ,z 1 ), 
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where Cp = 1/2 (iV c — 1/N C ) — > N c /2 in the large iV c limit, the factor 1/tt comes from the 
longitudinal integration measure, and 

K(x 2 0i x 21 ) = (1/4) [e(x 20 ) - <x 21 )} 2 . (2.4) 

However, unlike the situation of 3+1 dimensions, the kernel K{x 2 q,X2i) vanishes for x 2 
outside of the interval (cco,cci), (we assume that x\ > xq). This special property of one 
dimensional propagator implies that, due to color screening, no gluon can be emitted 
outside of the "parent" dipole in the transverse impact space. This fact is ultimately 
responsible for our final result that the total cross section does not grow with energy. 
Similarly, the probability of finding two soft gluons in the meson is 

</> (2) (x;, Zi) = \(K(x 32 , x 30 ) + K(xsi, x 32 ))\K(x 20 , £2l)</> (0) (£ , Xl ' Zl )' ( 2 - 5 ) 

where A = g 2 N c /(2ir). Since in the large- N c limit the color structure of a gluon can be 
identified with that of a quark-antiquark pair, ( |2.5| ) has a simple interpretation in a dipole 
picture: The second gluon at X3 is either emitted from the dipole consisting of one valence 
quark and the anti-quark part in the first gluon, or emitted from the dipole consisting of 
the valence anti-quark and the quark part of the first gluon. We emphasize that one could 
speak of splitting a gluon into a pair of quark and anti-quark only in a large- iV c language. 
In this next section, a clearer discussion on the relation of this picture to the usual large- iV c 
planar graph analysis will be provided. 

As in ||, we next introduce a generating functional $(xo,xi, zi,u(x, z)) whose nth- 
order coefficient in an expansion in u(x, z) is the probability amplitude of finding n soft 
gluons in the meson. Clearly, $ factorizes into (p(°\xo,Xi,Zi)Z(xo,Xi, u(x,z)). One can 
of course obtain these coefficients iteratively. More directly, it follows formally from ( |2.3j ) 
and ( |2.5| ) that Z satisfies the following integral equation: 

f Xl f Zl dz 

Z(xq,Xi,Zi,u) = 1 + A / dx 2 u(x2,Z2)Z(x2,Xl,Z2,u)Z(xo,X2,Z2,u). (2.6) 

Jxq J z z 2 

The structure of this integral equation can also be seen more clearly by examining figure 
2. Note that a cut-off zq to the longitudinal integration over z 2 has been introduced. 






Fig.2 

However, ( p?6|) as it stands is inconsistent, as pointed out in || for the analogous 
equation in the case of 3+1 dimensions. If we define <j)(°\xo, Xi, Z\) as the inclusive 
probability of finding the valence quark and anti-quark in the meson, then production 
of soft gluons should not change the probability. Consequently one has 4>(°\xq, x\, Z\) = 
$>(xo,Xi,zi,u = 1), or, Z(xo,Xi,Zi,u = 1) = 1. However, this condition is not satisfied 
by ( |2.6| ). To rectify the situation, we notice that the "bare" wave function (f>^> used in 
( EOl) must be replaced by (f)^ A(xq, xi, Z\) where the factor A takes into account virtual 
corrections to the emission of one real gluon. A similar factor must be used in each dipole 
constructed from emitted gluons. Since Z{xq,x\,z\,{S) = A(xq, X\, z±), it follows that it 
should appear as the inhomogeneous term in (|2.6|). Introducing rapidities y = \n(zi/zo), 
y' = ln(z2/zo), it is convenient to treat the soft gluon emission as an evolution process. 
Normalizing A(xq, x\, zq) = 1, it can be interpreted as the survival probability of no soft 
gluon emission at y > 0. It follows that ( |2.6| ) should be further modified by inserting a 
factor A(xq, x±, zxZq/z^), i.e., the probability of no gluon emission as one evolves from y' 
to y. 

This survival probability function satisfies a first order differential equation dA/dy = 
— Xfdx2K(x2o 1 X2i)A. It follows that A = exp(— Xxioy). The integral equation (|2.6|) is 
thus modified to become 

r xi r y i 

Z(x ,x 1 ,y,u(x,y))=e- Xxioy + \ / dx 2 / e~ Xxio(y - y ^dy'uix.y') 

Jxo Jo (2.7) 

Z(x 2 , xi, y', u)Z(x , x 2 , y', u). 



Note that the desired normalization condition at u = 1 is satisfied by (|2.7|) . There is no 
ultraviolet divergence in the virtual corrections in 2+1 dimensions since QCD is super- 
renormalizable, unlike the situation for 3+1 dimensions. 



3. One Pomeron exchange 

3.1. Topological structure of a meson state 

Knowing the integral equation for the generating functional Z, it is straightforward to 
write down the integral equation for the dipole density. Before doing so, let us first provide 
an intuitive picture of such a density function from the point of view of large- N c expansion 
where the dominant topological structure is that of closed-string exchange, figure 3. The 
cylindrical surface in figure 3 is understood to have been populated with arbitrary number 
of gluons, provided that they are connected in a "planar" fashion in their color structure. 
This can best be represented if each gluon is labeled in a "double-line" representation, i.e., 
as far as color is concerned, a gluon can be thought of as a pair of quark and antiquark. In 
this fashion, the cylindrical surface is "tiled" by polygons whose edges are gluon lines. The 
color structure is such that the boundary of each polygon corresponds to a closed "quark 
loop" , formed by identifying a quark line from each gluon of the polygon edges. 




B 



A' 




B' 



Fig.3 

Figure 3 can also be obtained by joining of a closed-string propagator, i.e., the long 
neck, to wave functions at both ends of the graph. By construction, this closed-string 
propagator corresponds to a color-neutral exchange. (In what follows, it will be represented 
by the color-neutral component of two-gluon exchange. Here we shall concentrate on the 
topological structure of the wave-functions at each end.) Each wave function, before the 
closed-string exchange, topologically has the structure of a disk. In order to join to the 
closed-string, a "window", or more precisely, a "quark-line" boundary must be cut out. In 
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a light-cone wave function approach, this quark-loop represents the propagation of a color- 
dipole. The "distribution" of this quark-loop can then be understood as the desired dipole 
density. The quark-gluon diagram and the corresponding surface diagram are illustrated 
in figure 4. 



Fig.4 



3.2. The dipole density in a meson state 

Let us now return to our treatment of the light-cone wave function of a fast meson. Let 
n(xo, Xi, x, x, y) be the dipole density in a meson with valence quark at xq and anti-quark 
at x\. The dipole has a separation x and center of mass position x. The density is defined 
for dipoles with the lowest longitudinal momentum greater than exp(— y)p+, where p+ is 
the total longitudinal momentum of the meson. By keeping track of the center of mass, 
our treatment here represents an extension of that treated in [Rl and H . 



From graphs in figure 5, it is easy to see that the integral equation governing the 
dipole density is 

n(x , xi, x, x, y) =e~ Xxwy 5(\x\ - x 10 )5(x - x 10 ) 



Xl 



+ X I dx 2 I dy'e Xxi0 ^ v Hn(x , x 2 , x, x, y') + n(x 2 , x 1 ,x, x, y') 



(3.1) 
where X\q is the center of mass position l/2(xo +Xi) of valence quarks. The delta function 
in the first term represents the contribution from the valence quarks. Again we assume 

Xi > Xq. 





+ 




+ 




To solve Eq. ( pT[) , use 



Fig.5 



duj 



n(x Q ,xi,x,x,y) = / -— e yu; n u ,(xo, x x , x, x), 

* Z7TZ 



(3.2) 



where the contour of the integral runs parallel to the imaginary axis and to the right of all 
singularities in n w . The integral equation for n^ becomes 

(u + \xio)n u (xo, xi, x, x) =5(|x| — xiq)S(x — xio) 

(3.3) 



+ A / dx 2 [n UJ (xo, x 2 , x, x) + n w (x 2 , xi, x, x)]. 

J x 



The above equation can be further simplified by taking a Fourier transform with respect 
to the center of mass x, 

n w (xo, xi, x, q) = / dxe %qx n LJ (xQ,xi,x, x). 



One of our key results is an exact solution to (13] 



-iqx 



w n LJ {x ,x 1 ,x,q) 



_5{\x\ -xiq) 
uj + Axio 
2A 



+ 0(xio - |x| 



4A 2 



'{u + X\x\) 3 q 



Xio — x 
sin | q 



(uj + x\x\y 



Xio — x 
cos I q 



(3.4) 



This solution can be understood intuitively. The first term is due to the dipole formed by 
valence quarks, the parent dipole. The second term comes from dipoles formed by emitted 
gluons. In particular, the step function in the second term tells us that the size of induced 
dipoles cannot be greater than x\q. This is a direct consequence of the fact that no gluon is 
emitted outside of the parent dipole. (We stress that n w (xo, xi, x, q) is not to be confused 
with a dipole density at a given transverse momentum q. Indeed, it is not even positive 
definite.) 

Transforming n^ in ( |3.4|) back to the rapidity variable, one finds 



e 



-zqx 



10 n(x , xi, x, q, y) = S(\x\- x 10 )e Xx ^+ 



4X 2 y 2 . .xio — \x\ xio — |x| \ (3-5) 



Of I h -Mx\yf 4X2 y 2 ■ { X 10 -\x\ X 10 -\x\ \ 

8{x 10 -\x\)e i |y I sm( q)+2Xycos{ q) I . 



Surprisingly, the transform of ( |3.5| ) back into the center of mass coordinate takes on a very 
simple form: 

n(x , xi, x, x, y) = 5(\x\ - x 10 )S(x - x 10 )e~ Xxioy + 6(x 10 - |x|)e~ A|x ' y [A 2 |/ 2 

(e(xi - x - |x|/2) - e(x -x+ |x|/2)) + Ay (8(x x - x - |x|/2) + S(x -x+ |x|/2))]. 

(3.6) 

It is easy to see that this expression for the density is positive definite, as it should. 

3.3. Meson-meson non-forward scattering 

Consider meson A collides with meson B, in their center of mass system. Denote 
the wave function of meson A by ^^(xio, z) exp(ipxi ) and the wave function of meson 
B by ^f b(x' 10 , z') exp(—ipx' 10 ). After integrating out their center of mass positions, the 
scattering amplitude between these two mesons can be expressed as 

A(s,t) =is dx 10 dx' 10 $ a(xio)$ b(x' 10 )A(x 10 , x' 10 ,Y, q) , (3.7) 

where $a(xi ) = / dz\^ A (x 10 , z)| 2 , $b(x' 10 ) = f dz'\^ B (x' 10 , z')\ 2 , and A(xx ,x / 10 ,Y,q) 
is a scattering amplitude between two mesons with fixed dipole sizes. Under a two gluon 
exchange approximation, it becomes 

/dk 1 f 

o~ u2( _ u)2 / dxdx'n(x ,x 1 ,x,q,Y/2)n(x' ,x' 1 ,x',-q,Y/2) 

sm(kx/2) sin((g — k)x/2) sin(/cx'/2) sin((g — /c)x'/2)), 

(3.8) 



where Y = Ins, and t = —q 2 - Thus, the amplitude given in Eq. (|3.7|) can be associated 
with a (hard) Pomeron exchange. Indeed as we shall see shortly, the leading contribution 
to A(s,t) in Eq. ( |3.7| ) conies from a logarithmic branch point in the complex angular 
momentum plane. 

Let us begin by first evaluating the following integral: 



/ 



dxn(xo, Xi,x, q, Y/2) sm(kx/2) sm((q — k)x/2). 



We first point out that this quantity depends on the center of mass position x\q only 
through the factor exp(zga?io) as one can easily seen from formula (|3.5|) . However, this 
dependence, together the similar factor arising from the second dipole density in (ft.8|) , will 
be removed when the integration over center of mass positions is carried out in arriving at 
( p.7|) . (The net result is a momentum conservation factor which has been dropped from our 
formula (|3.7| ).) We henceforth ignore the dependence on x\q in the dipole density. Next 
observe that, when Y large, the leading term in (|3.5| ) is proportional to Y 2 . It follows that 



/ 



dxn(xo, x\, x, q, Y/2) sin(kx/2) sin((q — k)x/2) = 2X 2 Y 2 

(3.9) 



%-fe)rv ( iax 10 /2 XY + iq 



{X 2 Y 2 + 2iXqY)(X 2 Y 2 + 2iXqY - Ak(q - k)) 



where 3(F) denotes the imaginary part. 

In arriving at (|3T9|), we have initially assumed X\ > xq. However, it is easy to see 

( |3.9|) depends only on |xio|- It follows that the corresponding integral in (|3.8|) involving 

the second dipole density can be obtained directly from ( |3.9|) by letting q — ► —q, k —* —k 

and x — > x' . Substituting these results into (|3.8|) and performing the integral over k, we 

find 

4 1 ,2Aycos(|xio|-|xi |)g/2) 



A(x w ,x' 10 ,Y) =- g 



q 2 X 2 Y 2 + Aq 2 



x2 v 2r ex P(^(l^lol + klol)/ 2 ) , „ r in 

1 (XY + 2iq) 2 (XY + iq) + JJ ' 
This amplitude is non-singular in the limit q — > 0. When \q\ « XY, it becomes 



(3.10) 



4g 4 
A(x 10 , x' 10 , Y) = —^ sin(^i /2) sin(ga; 10 /2), (3.11) 



and, in the limit q = 0, it behaves as 



A(x 10 ,x[ ,Y) = 9 * lx ^ X ' wl . (3.12) 
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We point out that the amplitude decreases as 1/Y at high energies. This result is 



in agreement with a calculation done in Ref. [p] , using an effective action approach. 
Note in the dipole picture, the important ingredient of virtual correction is derived rather 
indirectly. In [jTTf , however, the decrease of scattering amplitude is calculated with an 
effective action without additional input. ( |3.11| ) can also interpreted as composed of a 
Pomeron "propagator" 1/(XY) and its coupling to dipoles x%o and x' 10 . In this language, 
the coupling of a BFKL Pomeron to a dipole of size xio is given by 

P(q) = sin(gxio/2). 

Q 

Finally substituting Eq.( |3.10|) back to (|3.7|), one arrives at 



where 



$ A (q/2) = I dx$ A {x)e iqx/2 , $ B (g/2) = / dx<S> B {x)e iqx/2 
>o Jo 



and $yi(— q/2), $b(— q/2) are their complex conjugates. (We assume that both &a(x) 
and $s(x) are even functions.) For the forward limit q = 0, one can use (|3.12|) with [a?io| 



and |x' 10 | replaced by the effective sizes of meson A and meson B. Note Also that the 
high energy behavior ( |3.12| ) can be understood as due to a logarithmic branch point in the 
complex angular momentum plane, i.e., 

The branch point is located at uj = 0. 

Physically, the asymptotic decrease of the scattering amplitude is due to a suppression 
of the dipole density at large rapidity. This can be understood as follows. For 2+1 
dimensions, since there is only one transverse dimension, the emitted gluons in a meson 
is confined between the valence quarks, and the branching process of emission of soft 
gluons must compete with a suppression factor exp(—X\x\y) due to virtual corrections. In 
3+1 dimensions, this competition is won by the branching process, hence a BFKL Pomeron 
emerges to the right of J = 1 in the complex angular momentum plane. In 2+1 dimensions, 
however, the branching process losses the competition and the suppressing factor due to 
virtual corrections dominates, thus leading to a logarithmic branch point at J — 1. 
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4. Two Pomeron exchange 

4-1. Dipole pair density 

The next level contribution to the scattering amplitude in the leading logarithmic 
approximation is from exchanges of four gluons between two mesons. In order to calculate 
this contribution, we need to know the density of two dipoles in a given meson state. 
Let n 2 (xo, x\,x a , x a , Xf,, Xb, V) be the density of two dipoles with respective center of mass 
positions x a and x b and separations x a and Xb- Again this density is defined for dipoles with 
lowest longitudinal momentum greater than e~ y pj r , where p+ is the meson longitudinal 
momentum. (In deriving an integral equation for n 2 , we will first focus on its dependences 
on xq, x\ and y. In the following, for compactness, we shall sometimes drop its dependences 
on x a , x a , x b , and Xb- They will be re-inserted when necessary.) 

After a careful calculation, we have found that the leading term in the double density 
is simply the product of two single dipole densities. This will be demonstrated below. 
This surprising result implies that, to the leading order, there is no correlation between 
two single dipole densities. This is not the case for 3+1 dimensions. We take this as one 
of the remarkable features of the Pomeron physics in three dimensions. 

The integral equation can be written down readily by examining the graphical equation 

in figure 6: 

f Xl f y , 

n 2 (x ,x 1 ,y) -X dx 2 / e ~ Xxio{ - v ~ v 'dy'[n(x 2 , x 1 , x a ,x a , y')n(x , x 2 , x b , x b , y') 

Jx JO 

+ n(x 2 , xi, x b , x b , y')n(x , x 2 , x a , x a , y') + n 2 (x 2 , x\, y') + n 2 (x , x 2 , y% 

(4.1) 

where again we assumed x\ > xq. The first two terms come from finding two dipoles in 

separated dipoles, as indicated in figure 6, which serve as the inhomogeneous term in the 

integral equation. Note that the single dipole density depends on the separation of the 

dipole only through its absolute value, the same is true for the double dipole density, as 

evident from the above equation. 






+ \ " + 




Fig. 6 
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The single dipole density, as given in (3J5), is already lengthy. Fortunately, in order 



to calculate the leading term in the double dipole density, it is enough to only keep the 
leading term in the single dipole density as the input. This term, from (|3.6|), is 



(4.2) 



n(xo, xi, x, x, y) =X 2 y 2 e xy 9(xio — x)(e(xi — x — x/2) — e(xo — x + x/2)) 
= 2\ 2 y 2 e- Xxy 6(x 10 - x)9(x 1 -x- x/2)9(x - x - x/2). 



Plugging this into (|4.1|) , the integral equation reduces to 



n 2 (x ,£i,2/) -nf\x ,Xx,y) 

+ A / ^ dx 2 F ' dy'e-^ob-rt (n 2 (x 2 ,x 1 ,y') + n 2 (x ,x 2 ,y')) . 

Jx JO 



where 



n 2 



(0) 



> , x u y) = 4\ 5 e~ Xx ™y F y' 4 e x ^-^- x ^' dy'[(x a - x b - (x a + x b )/2) 

Jo 



9(x a - x b - (x a + x b )/2)9(xi - x a - x a /2)9(x b - xq - x b /2) + (a <-> b)}. 

(4.4) 

The presence of various step functions assures that both dipoles are found only within the 

parent dipole, and they do not overlap. 

Eq. (|4.3|) can be solved by taking a Laplace transform with respect to y, and this is 

carried in the Appendix. The leading order term in y takes on a surprisingly simple form 

n 2 (x , xi, x a , x a , x b , x b , y) =4X 4 y 4 e~ x{Xa+Xb)y 9(x a - x b - (x a + x b )/2) 

(4.5) 
9{x\ — x a — x a /2)9(x b — xq — x b /2) + (a <-> b). 



Comparing this result to the single dipole density in Eq. ( |4.2| ), we find that the leading 
term of the double dipole density, as given in ( |4.5|) , is just the product of two single dipole 
densities, the only additional requirement is that the two dipoles do not overlap. The first 
term corresponds to the situation that the center of mass position of dipole a is above 
the center of mass position of dipole b. The second term represents the other possibility. 
This result is not at all obvious before we solve the integral equation (|4.3|) . Note that the 
inhomogeneous term in ( |4.3|) , given in the first graph on the r.h.s. of figure 6, is different 
from ( |4.5| ), since its leading term in y is proportional to x a — x b — (x a + x b )/2, which is 
absent in ( f4.5|) . With this factor, the density is not homogeneous regarding to the center of 
mass positions. The result in ( f4.5| ) is rather homogeneous in the center of mass positions. 
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4-2. Meson-meson amplitude due to two Pomeron exchange 



With formula ( fO| ) of the dipole pair density, we can calculate the scattering amplitude 
due to four gluon exchanges between a pair of dipoles in one meson and another pair of 
dipoles in the other meson. The scattering amplitude with zero momentum transfer, for 
given parent dipole sizes xio and x' 10 , is given by 

2! f dk 
A 2 (x 10 x' 1Ql q = 0, Y) = ' / —dx a dx b dx' a dx b n 2 (x , xi, x a , k, x bl -k, Y/2) 

dk 
n 2 {x' , x[, x' a , k, x' b , -k, Y/2) [8g 4 a 2 d(k a , k, x a )d(k a , k, x' a )] (4.6) 

2%k^(k - k a y 

dkh 
[ 8 # 4 o ,2n i i \-> d(h, -k, x b )d{k b , -k, x' b )], 
2ivk^(k + k b ) 2 

where the term in the first bracket is from the exchange of two gluons between dipoles x a 
and x' a , and the term in the second bracket is from the exchange of two gluons between 
dipoles Xb and x' b . In ( |4.6| ), d(k a , k, x a ) = sm(k a x a /2) sm((k — k a )x a /2), and the dipole den- 
sity 77-2(2:0, x i5 x «! k, x b , —k, Y/2) is the Fourier transform of 712(2:0, x i> x ai x a , x b , x b , Y/2) 
with respect to the center of mass positions, (relative to the center of mass of the par- 
ent dipole). Finally, the numerical factor 2! is from possible pairing among two pairs of 
dipoles, and two factors 2! in the denominator are due to indistinquishability of two dipoles 
in a meson. A version of Eq.( f4.6j ) in two transverse dimensions has been derived in the 
Appendix of 0. 

Were the dipole density 712 more complicated than given in Eq. ( |4.5| ), it would be 



a hard task to perform the integral in (f4.6| ). Fortunately, Eq. ( }4.5|) tells us that in the 
leading order, the dipole pair density is the product of two single dipole densities, up to the 
requirement that the two dipoles do not overlap. Recall in our calculation of the amplitude 
due to one Pomeron exchange, the major contribution is from dipole size x ~ 1/(AY), which 
is very small at high energies. The smallness of dipole sizes then allows us to neglect the 
requirement that two dipoles do not overlap in a given meson, then the dipole pair density 
is simply the product of two single dipole densities. The integral in (|4.6|) factorizes. Thus, 
the two Pomeron exchange amplitude is given by 



1 f dk 

A2{xiQ,x' 1Q ,q = Q,Y)=- I —A(x 10 ,x' w ,k,Y)A(x w ,x' w ,-k,Y) 

2 J Zn 

= 2 / ^ A2 ( x w,x' 10 ,k,Y), 



(4.7) 
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where the last equality derives from the fact that the single Pomeron exchange amplitude, 
as given in (|3.10|) , is an even function of k. The above formula allows us the calculate the 
forward scattering amplitude A 2 from the knowledge of A. Performing the integral over k 
in (|4.7|) and keeping only the leading term, we have 



A 2 (xi ,a;i ,g = 0,Y) = ^ 2y2 I -\x 10 - x' 10 \ 3 + -(x 10 + x' 10 f - x 3 10 - (x' 10 fj . (4.8) 

So the amplitude is suppressed by the factor l/(lns) 2 , comparable to the sub-leading term 
in the single Pomeron exchange amplitude. The dependence of A 2 on the dipole sizes is 
cubic, as the correct dimensionality requires. In principle one can also compute A 2 for a 
finite momentum transfer. The formula is given by 

1 f dk 

A 2 (x 10 , x' 10 , q,Y) = - —A(x 10 , x' 10 , k, Y)A(x 10 , x' 10 , q-k,Y). (4.9) 



5. The triple Pomeron coupling 

We have found a simple formula for the scattering amplitude of two Pomeron exchange, 
indicating the simplicity of Pomeron physics in the three dimensional QCD. The next 
important object to calculate is the triple Pomeron coupling. Armed with knowledge of all 
these, one may attempt to construct the Regge field theory which helps one to study high 
energy scattering systematically. It is interesting, as one may expect, that once again the 
triple Pomeron coupling can be calculated exactly to the leading order, unlike the situation 
in the four dimensional QCD. 

The physical process involving triple Pomeron coupling is A + B^C + X, where A, 
B and C are mesons and X is anything. The generalized optical theorem states that the 
partial cross section of this so-called one-particle inclusive process is given by the imaginary 
part of amplitude of the elastic process A + B + C — > A + B + C [|12| . To be more precise, 
let p be the out-going momentum of C, then Ecda/d 2 p = 2s cL-'i 2 ^(^ 6 )' w h ere ^6 is the 
connected part of the forward scattering amplitude oiA + B + C^A + B + C. To extract 
information concerning the triple Pomeron coupling, we consider the case C = A and the 
one particle inclusive process is dominated by a Pomeron exchange. Instead of considering 
some definite meson states, let us consider mesons with a fixed dipole size. Let xio the 
dipole size of the B meson and x' 1Q the dipole size of the A meson. Furthermore, we define 
M 2 be the mass squared of X, rapidity y = ln(s/M 2 ), where s is the total energy squared. 
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Also define Y = ln(s/m T ), wit is the transverse mass of the outgoing A meson. With 
these notations, it is easy to see that 

_ da da 

Ec-pr- = s - 



d 2 p dqdM 2 ' 

where q is the momentum transfer. Consider large y, the process is then dominated by a 
single Pomeron exchange, so 






2 



where we have used the Pomeron "propagator" and its coupling to meson A, (3bx denotes 
the Pomeron coupling to B and X. Again applying the optical theorem, \J2x Pbx\ 2 
can be written as a product of the triple coupling coupling, the Pomeron propagator at 
Y — y = ln(M 2 /m T ) and the Pomeron coupling to meson B. Note that this coupling to 
B is the one at zero momentum transfer. Putting everything together, one has 

da M 2 8g 6 x 10 , . .. 2 .. , . 

Ec d^p ~ — \*y*(Y-y)\q\ ^^0/2)) V 3P ( q ), (5.1) 

where Vzp shall be referred to as the "triple Pomeron coupling" . 

Next, we calculate the 6-point amplitude Aq and use the optical theorem to relate 
this amplitude to the above quantity, and finally read off V^p. One can imagine that such 
amplitude is due to one Pomeron exchange between A and B and one Pomeron exchange 
between another A and B. So for meson B we need a dipole pair density. The rapidity 
available for one Pomeron exchange is y. However for the dipole pair in meson B, the 
rapidity available is different from y, and y serves only as a rapidity cut. We then follow 
to define the dipole pair density depending on both Y and y, here Y is the maximum 
rapidity. 

The desired equation again can be readily read off from a graphical equation, which 
we spare here: 

n 2 (x , xi, y, y) = A / dx 2 e~ Xxio( - y ~ y) [n(x 2 , xt,x a , x a , y)n(x , x 2 , x b , x b , y) 

fXi ry / ( 5 - 2 ) 

+ (a^b)} + \ dx 2 dy'e- Xxioiy - y \n 2 (x 2 ,x ll y' 1 y)+n 2 (x ,x 2l y' 1 y)}, 

Jx Jy 
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where all quantities with subscript a are those associated to meson A, and those with 
subscript b are associated to meson B. The solution to this equation can be obtained in a 
similar way as we solved Eq.( f4.1| ). We skip the details and directly write down the answer 



n 2 (x , xi, y, y) =4A Ve- A( " i)( ^ } - A(Xa+Xb)( ^ )/2 (za -x b - (x a + x b )/2) 

9(x a -Xb~ (x a + x b )/2)0(xi - x a - x a /2)9(x b - x - x b /2) + (a <-»• b). 

(5.3) 

We remark that this solution is exact with the input of the single dipole density given 

in ( |4.2| ) . Next we need to transform this density into the momentum space, which is the 

quantity needed in calculate the 6-point amplitude. Assuming q = q a = —q b and ignoring 

the dipole sizes in various theta functions, we find 



n 2 (x ,x 1 ,q,-q,y,y)= dx a dx b e iq( - Xa Xh ' ) n 2 (x 0l x ll x ai x bl y 1 y) 
= 4\ 5 y 4 e- Xix « +Xb)iy+ y )/2 \— ^ ( e -(A( I /-S)+ig)x 10 + x) (5 4) 

(Ky-y) + t( i) ' 

+ 777 1 (e- (x(y -y )+iq)xw - 1) + c.c.]. 

(Ky - y) + w 

This formula is still a bit too long. Indeed the physical situation corresponds to X(y — 
y)xio >> 1, and this simplifies the above formula a lot: 

n 2 (xo, x t , q, -q, y, y) = 8A 5 ^io 7^ ~^ lf 2)2 e- X{Xa+Xb){y+ " )/2 - (5-5) 

Now we are ready to calculate the 6-point amplitude. By definition, this amplitude is 
given by 

dx a dx b dx' a dx' b n 2 (x , x u x a , q, x h , -q,Y - -, -)n(x' w , x' a , q, -) 

n(x[ , x' b , q, -)[8g 4 - . 2 d{k a , k, x a )d(k a , k, x' a )] (5.6) 

z Z7iK a {q — K a ) 



dk b 
2nkl(q~- h) 2 ' 



[8^ 4 2 . T-^d{kh, -k, x b )d{k b , -k, x' b )], 



where Y — y = In M 2 and other notations are the same as in the previous section. Recall 
that the single dipole density in the momentum space is 



n(x' 10 , x' a , q, y) = A\ 2 y 2 -sin(qx' 10 /2)e x <\ 
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again we ignored the dipole separation as compared to x' 10 . Substituting this formula and 
the formula in Eq.(|575|) into (pT6|), we find 

A 6 (q) = -A M y (A2(y _- )2 + g2)2 ^r (sm(^ 10 /2)) 
(V J dx a dx' a e- x *" Y / 2 - x <*/ 2 2nk J* a _ K)2 d(K, q, x a )d{k a , q, x' a )\ . 

Finally performing the integral in the above equation and one is led to a quite lengthy 

formula 

MO) =2*9 8 M 2 \ 7 f - * 10 (sm(qx' 10 /2)) 2 

X 2 (Y-y) 2 -q 2 (5 - 7) 



(X 2 (Y - y) 2 + q 2 ) 2 (X 2 Y 2 + q 2 ) 2 (\ 2 y 2 + q 2 ) 2 ' 

This formula is actually not valid when q ~ Ay, since we have ignored the dipole size in 
the sine factors in the single dipole density as well as in the dipole pair density. This is 
to be contrasted with the case for one Pomeron exchange amplitude, (|3.10| ), where the 
dependence on dipole sizes is not ignored. (If one so desires, one can of course take the 
dependence on dipole sizes into account, leading to a much more complicated function of 
q.) When |g| « XY where Eq. (|5.7|) is valid, it simplifies to 



Md) = 2 ^ M \3 YHY + l]°2 {Y _y ) 2 ^ (M**io/2)) a • (5-8) 

Upon using the optical theorem, E c da/d 2 p = 2s( ^ 7i .p S ; (A 6 ), and comparing ( |5.8|) to (|5.1|) , 
we find that the triple Pomeron coupling depends on rapidities: 

a 2 v 6 
Vo P ~ 2 z . (5.9) 

3P \q\YHY + y) 2 (Y-y) l ] 



It tends to zero when both Y and y are large. It is also singular when y—t = \q\ — >• 0, just 
as in the 3+1 dimensional case |§. However the origin of this singularity in the present 
case is purely kinematic, as Eq. (|5.1[) shows. If the one particle inclusive cross section 
factorizes, one expects that V%p is independent of rapidities. Our result indicates that the 
inclusive cross section does not factorize and the triple Pomeron coupling vanishes as7 _1 , 
with Y — y fixed. 
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6. Discussion 

As we have seen, one remarkable property of high energy scattering processes in our 
toy world, 2+1 dimensional QCD, is the fact that no soft on-shell gluon can escape from 
hadron emitting it. This is due to the quantum interference of quark and anti-quark. This 
property is also the reason underlying the fact that the total cross section decreases at high 
energies, as was first pointed out in ||11|| . It is also responsible for the interesting fact that, 



in the leading order, there is no correlation between a pair of dipoles in a meson state, so 
the amplitude due to two Pomeron exchange is simply given by Eq. (|4.9|) . Physically, one 
expects that this holds true even for multi-dipole density. Thus, in the leading order, the 
amplitude due to n-Pomeron exchange is given by 

1 f n ~ X ( dk- \ 

A n (x 10 , x' w , q) = — Y[[ -j 1 A(x 10 , x'icn h) J A(xiq, x' 10 , q-^h), (6.1) 

i=l i 

where A(xio,x'i , k) is the amplitude due one Pomeron exchange between dipole x±o and 
dipole x' 10 . ( |6.1] ) is a straightforward generalization of (|4.9| ). Going to the impact parameter 



space, Eq. ( |6.1| ) implies that the eikonal formula for summing over all possible number of 
Pomeron exchanges is valid in the leading order. However, as can be seen in Sec. 3, there 
are sub-leading terms on each level of Pomeron exchanges, and these terms are comparable 
to the leading term on the next level of Pomeron exchanges. 

Of course, to incorporate the whole complexity of high energy scattering, one needs 
to include various couplings among Pomerons. The triple Pomeron coupling is already 
calculated in this work. It is interesting to see whether one can develop a Reggeon field 
theory based on our results. 
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Appendix 

Eq. ( |4.3|) can be simplified by a Laplace transform. Denoting the transform for n 2 by 
^2,w(^o, xi), where the variable u plays the role of a "complex angular momentum ", the 
integral over y' can be undone. Let 

n 2)CU (x , xi) = - — -j (u + X(x a + x h )) n 2 ,u(xo, xi), 

we obtain the integral equation for n 2 ^: 

(a; + Aaio)n2,w(#o, x\) = A 5 [(x a - x b - (x a + x b )/2)9(x a - x b - (x a + x b )/2) 

9(x\ - x a - x a /2)9(x b - xq — x b /2) + (a <->• b)] + X / dx 2 (n 2 ^(x 2 , x\) + n 2jLU (x , xi)) . 

Jxo 

(A.l) 
To solve (A.l), we use the following ansatz 

fi2,Lu(xo, xi) =X 5 (x a - Xb - (x a + x b )/2)9(x a - x b - (x a + x b )/2)9{x 1 - x a - x a /2) 

9(x b -Xq- Xb/2)fu(XQ, Xl, X a , X a , Xb, Xb) + (d<-+b). 

(A.2) 
The unknown function f u then satisfies 



rx b —x b /2 «i 

(u + Xxio)fu(xo,xi) = l + X[ dx 2 f UJ (x 2 ,xi) + dx 2 f U} (x G ,x 2 )\. 

J xn J x n +x n /2 



i-x b —x b /2 Mi 

dx 2 fu;(x2,Xl) + / 
'%0 Jx a +x a /2 

The solution is 



fcj = [V + X(x a ~Xb + (X a + X b )/2)] 1 , 

which is independent of xq and X\. Substituting this result in (A.2) and recalling the 
definition for n%^, we find 

n 2 ^(x , xi) = 4 x 4\X 5 [uj + X(x a + x b )]~ 5 [u + X(x a - x b + (x a + x b )/2)}~ 1 (x a - x b 

- (x a + x b )/2)9(x a -Xb- (x a + x b )/2)9(xi - x a - x a /2)9(x b - x - Xb/2) + (a <->• b). 

(A3) 

It is easy to transform n 2)UJ back to arrive at n 2 (xo,xi,y). There are many terms each 

proportional to a certain power of y. The leading term takes on the following simple form: 

n 2 (x Q , xi, x a , x a , Xb, Xb, y) =4X 4 y 4 e~ x( - Xa+Xb)y 9(x a - x b - (x a + x b )/2) 

(AA) 
9(xi — x a — x a /2)9(x b — xq — Xb/2) + (a <-»■ b). 
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Figure captions 

Fig.l: A soft gluon emitted either from the quark or the anti-quark of a dipole. The dashed 

line denotes the gluon. 

Fig. 2: A graphical illustration of Eq. (2.6) for the generating functional Z. 

Fig. 3: A string picture of meson-meson scattering by a closed string, a Pomeron, exchange. 

Fig. 4: Quark-gluon diagram of the process of two gluon exchange and its corresponding 

surface picture, obtained by replacing a gluon line by two quark lines. 

Fig. 5: A graphical equation for the dipole density, again dashed lines are soft gluon lines. 

Fig. 6: A graphical equation for the dipole pair density. 
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